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We present the results of a search for dark matter production in the monojet signature. We analyze 
a sample of Tevatron pp collisions at TeV corresponding to an integrated luminosity of 

6.7 fb _1 recorded by the CDF II detector. In events with large missing transverse energy and 
one energetic jet, we find good agreement between the standard model prediction and the observed 
data. We set 90% confidence level upper limits on the dark matter production rate. The limits are 
translated into bounds on nucleon-dark matter scattering rates which are competitive with current 
direct detection bounds on spin-independent interaction below a dark matter candidate mass of 
5 GeV/c 2 , and on spin-dependent interactions up to masses of 200 GeV/c . 
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The cosmological abundance of dark matter (DM) is 
now precisely known through the observation of its grav- 
itational interactions Yet the nature of DM itself 
remains a mystery, with many models of physics be- 
yond the standard model (SM) proposing DM candi- 
dates. Perhaps the best motivated DM candidate is a 
new weakly interacting massive particle (WIMP) with 
mass of 0(1 - 1000) GcV/c 2 . This class of DM candi- 
dates appears in many models of new physics with in- 
teractions that allow for DM detection through WIMP- 
nucleon scattering in direct detection experiments Q ■ 

While there is no conclusive evidence for WIMP- 
nucleon scattering, several recent direct detection ex- 
periments have yielded results suggestive of a low-mass 
(-10 GeV/c 2 ) WIMP [$% In light of these results, 
there has been significant interest [fj|-|9| in the potential 
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of collider searches to cither observe the production of 
DM particles (x), or to constrain the DM production 
rate. The collider mode of production that is expected 
to yield the most stringent bounds on the DM production 
rate is the monojet (pp — > XX+jet) mode, where the jet 
has a transverse energy of 0(100) GeV and originates in 
initial state radiation. Previous studies of the monojet 
signature have been performed (Iol - [l2l | in the context of 
searches for large extra dimensions. 

In this Letter, we present the results of the first di- 
rect search for collider production of DM in the monojet 
mode. We consider several models of DM production 
that are relevant for direct detection experiments. We 
assume x is a Dirac fermion [HI], and that production is 
mediated by a massive state which couples to DM and 
SM quarks. For this analysis, wc consider three models 
of dark matter production which consist of vector (Oy), 
axial- vector (Oav)> an d i-channel operators (O t ) as de- 
fined in 0, 0). A universal sum over all quark flavors 
is assumed for these operators. In direct detection ex- 
periments, the vector operator leads to spin-independent 
DM scattering, while the axial vector operator is spin- 
dependent. The ^-channel operator includes both spin- 
dependent and spin- independent terms. By considering 
the three types of operators, we are able to constrain 
both spin-independent and spin-dependent DM-nucleon 
scattering cross sections. 

In direct detection experiments, scattering rates are 
described by an effective theory containing DM in ad- 
dition to SM fields. As the momentum transfer in DM 
scattering is far lower than the mass of the particle medi- 
ating the interaction, an effective theory provides a valid 
description. In a collider environment, with large mo- 
mentum transfers, the effective theory approach is not 
necessarily valid and may change the predicted cross- 
section and kinematics of the DM model [1, 0, E3] ■ We 
thus consider two possibilities: (1) that these contact in- 
teractions are also a good description of collider DM pro- 
duction, and (2) that the production of DM at the Teva- 
tron proceeds through the exchange of a new particle. 
The new mediator particles which lead to the operators 
Oy, O av , and Ot, are a heavy vector, axial- vector, and 
a scalar "squark" respectively. When constraining case 1 
we implement models of DM production with very heavy 
mediators, well above the Tevatron reach (at 10 TeV), 
while for case 2 we consider light mediators, within the 
kinematic reach of the Tevatron (100 GeV/c 2 for Oy and 
Oav and 400 GeV/c 2 for O t ). 

Wc perform the analysis utilizing 6.7 fb _1 of Tevatron 
pp collisions at y/s=1.96 TeV recorded by the CDF II 
detector. The CDF II detector is described in detail else- 
where [lij and consists of tracking systems immersed in 
a 1.4 T magnetic field, surrounded by calorimetery that 
provides coverage for \r)\ < 3.6 [lit ]. A system of drift 
chambers external to the calorimetery provides muon de- 
tection capability for \rj\ < 1.5. 



The DM candidate is expected to interact minimally 
with the CDF II detector resulting in large missing trans- 
verse energy {$t) [13 • We analyze a sample of events 
consistent with this characteristic, as collected by a $t 
online event-selection (trigger) algorithm which selects 
collision events with $t > 40 GeV. We find that the 
trigger has a 90% selection efficiency for events with J^t = 
60 GeV, rising to an efficiency of 95% for E T > 70 GeV. 

The selected events are required to have been recorded 
with fully-functioning calorimeter, muon, and tracking 
systems. We require events to have $t > 60 GeV. Wc 
reconstruct jets from calorimeter energy deposits using 
a cone algorithm fl8j with a radius in pseudorapidity- 
azimuth space of 0.4. The jet energies are corrected for 
variations in detector response and instrumentation, and 
the extra contribution from additional pp pair interac- 
tions in the same event [l9| . Events are required to have 
exactly one jet with Et > 60 GeV and |?7| < 1.1. To re- 
ject events arising from non-collision sources, we require 
significant track activity within the jet cone. Events in 
which the jet does not contain at least one track with 
a transverse momentum {pt) of at least 10 GeV/c arc 
rejected. We reject events in which the jets are re- 
constructed in a partially instrumented regions of the 
calorimeter. To remove photons we require that the elec- 
tromagnetic fraction of the total energy deposited in the 
calorimeter systems to be below 0.85. Similarly, to re- 
move events with muon bremsstrahlung from cosmic rays 
or beam-detector interactions [io| we require an electro- 
magnetic fraction of greater than 0.35. To accommodate 
extra jets arising from initial state radiation, we retain 
events with one additional jet with an Et of less than 
30 GeV and \rj\ < 2.4. 

The sample of events that pass the above selections 
is dominated by background contributions from QCD 
multijet processes in which one (or more) of the jets is 
mis-reconstructed. Improper reconstruction of a jet pro- 
duces an event topology in which the J^t is aligned with 
the mis-reconstructed jet. To reduce this background 
wc require a minimum separation in azimuthal angle of 
A(/)(^t, jet) > 0.4 between the direction of J^t and that 
of any jet with Et > 20 GeV in the event. We also 
require a separation of A<fi($T, jet) > 2.5 between the 
direction of J£t and the leading jet. We achieve further 
reduction of the multijet background to our search by uti- 
lizing an artificial neural network (NN) designed to sepa- 
rate multijet events from electroweak processes. The NN 
combines event quantities including the separation in az- 
imuthal angle between jets and $t, jet energies, $t, and 
the number of jets, returning a single numerical value for 
each event. In training, the NN was optimized to isolate 
simulated Z and W boson events from a sample of data 
events in which the most energetic jet had Et < 60 GeV, 
or in which there were more than three jets. We find that 
approximately 85% of multijet events produce a NN value 
of less than 0.3, and reject these events. 
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In the remaining sample of events, we expect signif- 
icant contributions from Z and W boson processes, in 
which the Z or W decays leptonically. We reduce the con- 
tribution from these processes by vetoing events which 
contain one or more tracks with pr > 10 GeV/c that are 
not embedded within a jet. 

The events passing the above selections form our anal- 
ysis sample, and are examined for the presence of events 
arising from DM production. Within this sample we ex- 
pect significant contributions from Z boson processes in 
which the Z boson decays invisibly to neutrinos. In ad- 
dition, we expect W boson processes to contribute when- 
ever the lepton from the leptonically decaying W boson 
is outside of the acceptance of the CDF tracking system. 
We model W and Z boson contributions to our anal ysis 
sample using simulated events generated by Alpgen 21 1 
with Pythia [22| for particle showering and hadroniza- 
tion. The Z and W boson background contributions are 
determined assuming NNLO calculations [23| of the in- 
clusive production rates. 

Minor backgrounds include tt modeled with Pythia, 
and single-top processes modeled with MadGraph [24[ 
plus Pythia for particle showering and hadronization. 
A top-quark mass of 172.5 GeV/c 2 is assumed for the 



tt [25|, l26| s-channel, and f-channel [27j processes with 
cross sections of cross sections of 7.04, 1.05, and 2.1 pb, 
respectively. We account for the expected diboson 
(WW,WZ,ZZ) [HI contributions to our selected sam- 
ple with a Pythia simulation, and normalize the rates 
of the WW, WZ, and ZZ processes to 11.34, 3.47 and 
3.62 pb. All simulated samples in this anal ysis include 
a detailed GEANT-based detector simulation [291 ] and as- 
sume CTEQ5L [30| parton distribution functions. 

While our NN requirement rejects the main multijet 
contamination, we model the remaining multijet back- 
ground using reweighted data events. We determine the 
likelihood of an event in our analysis sample to have orig- 
inated from a multijet background process by utilizing a 
sample of events with relaxed kinematic selections such 
that events with any number of jets with Et greater than 
35 GeV are accepted. Events meeting this relaxed selec- 
tion constitute the derivation sample. To maintain exclu- 
sivity, we remove all events entering the analysis sample 
from the derivation sample. In addition, we denote the 
probability that a given event originated in a multijet 
process as the multijet probability (MJP). 

The MJP is obtained as the fraction of events in the 
derivation sample remaining after subtraction of all sim- 
ulated backgrounds. We parameterize the MJP using six 
observables in order to mimic the topological characteris- 
tics of multijet events. These are the number of jets, 
the minimum separation in azimuthal angle between T£t 
and a jet, the ratio of the scalar sum of jet Et to its 
sum with the $t, the magnitude of the momentum im- 
balance from tracks with pt > 10 GeV/c, and the $t 
significance |3lj |. 



A given event in the analysis sample is assigned a 
weight by the MJP, as determined by its values of the 
six observables. The multijet background is modeled as 
the weighted sum of all events contributing to the anal- 
ysis sample. We find that the above method accurately 
determines the shape of the multijet background in all 
observables of interest. To obtain an appropriate nor- 
malization for the multijet contribution, we require that 
the sum of the number of events predicted by simulation 
and by the multijet prediction equal the number of data 
events observed in the sideband region which is defined 
such that the NN value is between 0.2 and 0.3. 

To test the performance of our data model, we form 
two additional samples that are exclusive of the analysis 
sample. We define an electroweak sample that is com- 
posed of events that pass all analysis selection criteria 
but have one or more tracks with pt > 10 GeV/c, that 
are not embedded within a jet. In addition, we define a 
multijet sample of events passing all the analysis selec- 
tion requirements except that they have a NN value less 
than 0.3, $t aligned with a jet, or have more than 2 jets. 
We find good agreement between the data and the SM 
prediction in both control samples. The Et distribution 
of the leading jet is displayed in Fig. [T] 

We model the potential contribution to our analysis 
sample from a DM signal of pp — > xx+jet using a Mad- 
Graph [HI generator that is interfaced with Pythia 
for showering and hadronization. We generate variants 
of the signal models, discussed previously, assuming DM 
masses between 1 and 300 GeV/c 2 . We find an efficiency 
of approximately 2% when imposing the analysis sample 
on simulated DM events. 

Systematic uncertainties affecting the normalization of 
simulated background components and DM signal arise 
due to uncertainty in the integrated luminosity (~6%), 
measured jet energy scale (~7%), parton distribution 
function uncertainties (~2%), efficiency of the trigger 
used for data collection (~2%), choice of the renormal- 
ization scale (~2%), and the amount of initial and final 
state radiation (~1%). In addition, a 50% uncertainty 
is placed on the normalization of the multi-jet predic- 
tion. The normalization uncertainties for the top 25l- 2l\ 
, Z/W [H], and diboson [H] processes are 10%, 8%, 
and 6% respectively. We include the effect on the jet 
energy scale uncertainty on the shape of observed quan- 
tities, and find this to be the dominant uncertainty. We 
include an uncertainty on the shape of the multijet back- 
ground, based on the observed variation in the multijet 
prediction between the analysis and electroweak samples. 

The total numbers of observed and expected events 
in the control and analysis samples is listed in Table |U 
In the analysis sample, we observe 52633 events which 
agrees well with the expectation of 53906 ± 6022. As we 
do not observe a significant excess over the number of 
events predicted by our background model, we proceed 
to quantify the maximum allowed DM production cross 
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FIG. 1: Distribution of the jet Et for the multijet control (a), electroweak control (b), and analysis (c) samples. The last bin 
contains the overflow. For the analysis sample, the jet Et of a representative signal process (xX+j e t) is shown normalized to 
the 90% confidence level upper limit production rate of 5.9 pb. 



section. 

We set limits on the DM production rate using a 
Bayesian likelihood [33| formed as a product of likeli- 
hoods over bins in the analysis region of the jet Et dis- 
tribution. We assume a flat prior on the signal rate, and 
a Gaussian prior for each systematic uncertainty includ- 
ing those affecting sample normalizations and shapes. 
We set Bayesian 90% confidence level upper limits on 
(j(pp — > XX + jet) f° r each of the models considered. The 
expected upper limits at each model point arc derived by 
randomly generating a series of pseudo-datasets, derived 
from the background prediction, and computing the me- 
dian of the distribution of resulting upper limits. The 
upper limits are listed in Table [TTJ We proceed to con- 
vert the limits into constraints on the DM-nucleon cross 
section following @, @| . A comparison of the CDF lim- 
its to several direct detection results is shown in Fig. [2] 
The CDF limits assuming light mediators are also shown. 
The CDF bounds extend beyond the experimental reach 
of direct detection searches, which are insensitive to DM 
with a mass of approximately 1 GcV/c 2 . For a DM mass 
of 5 GeV/c 2 , CDF bounds on spin- independent inter- 
actions are O(10 -38 ) cm 2 and are similar to the limits 
reported by the DAMIC [3^ collaboration. In the case of 
spin-dependent interactions, we report stronger bounds 
of O(10 -40 ) cm 2 for a DM mass of 1 GeV/c 2 , rising to 
0(1CT 39 ) cm 2 for a mass of 200 GeV/c 2 . 

In conclusion, we have performed the first collider 
search for DM in the monojet production mode. We 
have set limits on the DM production rate, and have 
constrained the spin-independent nucleon-DM scattering 
rate for a DM mass of roughly 1 GeV/c 2 , and between 1 
and 200 GeV/c 2 for spin-dependent interactions. 
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TABLE I: Event totals in the analysis and control samples. 
Uncertainties are systematic only. 



Source 


Multijet 


Electroweak 


Analysis 


Z 


6949 ± 840 


1280 ± 155 


22191 ± 2681 


W 


14986 ± 2007 


5582 ± 747 


27892 ± 3735 


Multijet 


165479 ± 82740 


1066 ± 533 


3278 ± 1639 


Other 


2194 ± 233.4 


149 ± 10.7 


545 ± 39.3 


Total model 


189608 ± 82787 


8076 ± 1011 


53906 ± 6022 


Data 


188361 


7942 


52633 



TABLE II: Expected (Exp.) and observed (Obs.) 90% C.L. 
upper limits (in pb) on the cross section of pp — > XX + j e t 
for the three operators (defined in text) Oav, Ov, and Ot, 
assuming contact interactions. The ±lcr variations on the 
expected limits are also shown. 





-Oav - 


-Ov - 




o t - 


M x 


Obs. Exp. 


Obs. Exp. 


Obs. 


Exp. 


(GeV/c 2 ) 


(pb) (pb) 


(pb) (pb) 


(Pb) 


(Pb) 



1 


5.9 


7 6+ 4 - 9 

' - D -3.6 


9.1 




6.5 


7.2 


5 


6.9 




5.4 




7.5 


8.1 


10 


4.5 




6.4 




5.0 


7.0 


50 


3.4 


7 0+ 4 - 2 

1 - u -3.4 


6.6 


7 c:+4.6 
' -°-3.4 


5.3 


6.5 


100 


4.5 


6-oi 3 2 6 9 


5.7 


6 2+ 37 


4.6 


6.1 


200 


4.8 


5.6± 3 2 2 r 


3.9 


5.6+1® 


4.2 


4.8 


300 


3.1 


6 1 +3 - 9 

0.±_ 2 7 


3.5 


5-6+ 3 2 1 


2.7 


5.1 



+4.3 
-3.6 
+4.4 
-3.6 
+4.2 
-3.1 
+4.0 
-3.1 
+3.7 
-2.9 
+3.2 
-2.2 
+3.2 
-2.4 
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FIG. 2: Comparison of CDF results to recent results from DAMIC 0, CoGeNT Q, XENON-100 [H, SIMPLE and 
COUPP [33| • Spin- independent (left) and spin-dependent (right) bounds are shown for the operators (defined in text) Oav, 
Ov, and Ot, assuming contact interactions. For comparison we also display CDF bounds assuming light mediators. 
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